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ABSTRACT 

Sublimation  from  blowing  snow  is  estimated  from  the  balance 
of  heat  and  mass  transfer  on  a  volume  of  air  and  nonuniform  ice 
particles.  Size  and  the  distribution  of  sizes  are  important  particle 
factors,  but  fragmented  and  abraided  shapes  appear  to  increase 
sublimation  only  slightly  above  the  estimate  for  spherical  parti- 
cles. Humidity  and  temperature  are  the  overriding  environmental 
factors,  but  atmospheric  pressure  and  solar  radiation  are  also  im- 
portant. The  effect  of  turbulence  on  convection  around  the  par- 
ticles does  not  appear  important  in  view  of  the  small  particle 
sizes  encountered.  Rather,  turbulent  transfer  determines  the  snow 
concentration  profile,  and  the  gradients  of  heat  and  water  vapor 
necessary  to  balance  the  sublimation  process. 

Keywords:   Snow  density,  snow  samplers,  weather  control,  evap- 
oration, sublimation. 
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Sublimation  of  Wind  Transported  Snow  -  A  Model 


R.  A.  Schmidt,  Jr. 


Although  a  small  process  in  the  total  hydro- 
logic  cycle,  sublimation  of  solid  precipitation 
during  redistribution  by  wind  may  be  a  sig- 
nificant part  of  the  water  balance  in  certain 
regions.  The  objective  of  this  work  is  to 
formulate  a  mathematical  model  to  estimate 
the  mass  of  snow  that  returns  to  vapor  as  the 
air-snow  mixture  moves  downwind  over  a  hori- 
zontal snow  surface.  2 

Following  a  brief  description  of  the  phenom- 
enon, sublimation  of  snow  is  considered  in 
three  parts.    Part  I  uses  an  existing  equation 
for  sublimation  of  an  ice  sphere  to  evaluate 
the  effects  of  meteorological  factors,  and  a 
modification  demonstrates  the  influence  of  par- 
ticle shape.    Part  II  shows  that  natural  snow 
concentrations  are  not  large  enough  to  cause 
significant  particle  interaction,  so  that  the  re- 
sults of  Part  I  apply  to  a  cloud  of  snow  par- 
ticles as  well.    Part  III  combines  this  knowl- 
edge with  expressions  for  the  vertical  gradients 
of  windspeed,  snow  concentration,  temperature, 
and  humidity,  to  estimate  sublimation  from  a 
vertical  column  above  a  horizontal  surface. 
Discussion  of  this  model  leads  to  several  inter- 
esting consequences  which  should  be  grounds 
for  fruitful  experimentation. 

The  title  employs  "sublimation"  to  describe 
the  process  under  discussion.  Although  sub- 
limation can  be  used  to  denote  either  deposition 
or  vaporization,  the  phase  change  from  solid 
to  vapor  is  of  primary  concern  here.  The 
process  is  assumed  to  be  quasi-stationary,  and 
sublimation  estimates  are  instantaneous  values 


My  initial  interest  in  this  problem  was  generated 
by  discussions  with  R.  D.  Tabler.  M.  D.  Hoover  and 
M.  Martinelli,  Jr.  helped  at  each  stage  of  the  develop- 
ment. J.  L.  Kovner  and  U.  Radok  reviewed  the  manu- 
script with  great  insight. 


that  do  not  take  into  account  the  change  in 
surface  area  of  the  particles  as  the  process 
proceeds.  Philip  (1964)  has  shown  that,  while 
not  fundamentally  correct,  such  analysis  pro- 
vides useful  approximations  for  many  meteoro- 
logical problems.  Calculations  are  included 
both  to  provide  more  meaning  to  the  equations 
and  to  given  values  that  help  keep  the  numer- 
ous factors  in  their  proper  perspective.  In  most 
cases,  these  results  are  presented  graphically. 
The  symbols  are  defined  as  they  are  intro- 
duced, and  are  listed  in  the  appendix 
for  reference. 

The  fact  that  the  snow  particles  vaporize, 
either  partly  or  sometimes  completely,  during 
transport  by  wind  is  not  directly  observable, 
since  measurement  problems  are  not  yet  solved, 
even  for  the  relatively  simple  situation  of  trans- 
port across  a  horizontal  surface.  At  windspeeds 
sufficient  to  move  snow,  natural  airflow  is 
turbulent,  and  windspeed  increases  with  height. 
Snowflakes  formed  during  precipitation  are 
quickly  reduced  toward  spherical  shapes  by 
mechanical  abrasion,  so  that  the  windblown 
particles  become  like  sand.  As  a  particle  leaves 
the  snow  surface  to  become  suspended  in  the 
turbulent  air  stream,  a  boundary  layer  develops 
around  it.  If  the  vapor  pressure  of  water  in 
the  surrounding  atmosphere  is  different  from 
the  saturation  value,  a  vapor  pressure  gradient 
forms  across  the  boundary  layer  and  sublima- 
tion begins.  The  rate  of  the  process  depends 
on  how  rapidly  air  is  exchanged  in  the  particle 
boundary  layer.  Energy,  in  the  form  of  heat, 
is  transferred  to  or  from  the  particle  surface 
by  conduction,  convection,  and  radiation.  Water 
vapor  is  transferred  by  diffusion  and  convection. 
These  same  mechanisms  operate  in  a  cloud  of 
snow  particles,  since  the  interparticle  spacing 
far  exceeds  the  depth  of  particle  boundary 
layers  under  natural  concentrations  of  snow 
suspended  in  turbulent  air. 
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PART  I 

SUBLIMATION  OF  A  SINGLE  PARTICLE 


A.  An  Ice  Sphere 

Two  laws  of  physics  provide  the  starting 
point  for  theoretical  considerations  of  the  sub- 
limation process.  Conservation  of  mass  provides 
a  relationship  in  the  form  of  the  diffusion 
equation.  The  relationship  for  heat  transfer 
comes  from  the  law  of  conservation  of  energy. 


1.  Diffusion  and  Conduction 

For  mass  transfer,  the  cla 
(Houghton  1933)  is  based  on 
equation  relating  the  rate  of 
(dm/dt)of  a  spherical  particle  ( 
surface  area  of  the  particle,  the 
of  water  vapor  in  air,  and  the 
difference  at  the  particle  surface 
equation  for  evaporation  of  a 
air  is 

(dm/dt)  =  4iTDr(p  -pr) 


ssical  approach 
the  differential 

loss  of  mass 
radius  r )  to  the 
diffusivity  (  D ) 
vapor  pressure 
.  The  resulting 

sphere  in  still 

[1] 


where  p  is  the  water  vapor  density  in  the 
remote  environment  and  pr  is  the  value  at 
the  particle  surface.  This  equation  expresses 
the  rate  of  sublimation  (or  condensation)  of  a 
spherical  particle  by  diffusion  without  ven- 
tilation. 

Heat  transfer  by  conduction  to  a  sublimating 
spherical  particle  depends  on  the  surface  area, 
the  temperature  difference,  and  the  thermal 
conductivity  of  air,  K.  The  expression  is 


Ls   (dm/dt)   =  4irKr(Tr  -  T) 


[21 


where    ls  is  the  latent  heat  of  sublimation. 

Temperature  in  the  remote  environment  is 

T,  and  t  is  the  surface  temperature  of  the 
particle. 

2.  Convection 

When  a  spherical  particle  sublimates  in  a 
steady  air  stream,  the  rate  is  higher  and  mass 
transfer  varies  around  the  particle  surface. 
Frossling  (1938)  predicted  this  pattern  from 
boundary  layer  theory,  and  verified  his  cal- 
culations with  measurements  of  mass  transfer 
around  a  naphthalene  bead.  The  rate  was 
maximum  on  the  exposed  face  of  the  sphere, 
and  decreased  to  a  minimum  at  the  point  where 
the  boundary  layer  separated  from  the  particle 
surface.  A  smaller  maximum  occurred  in  the 
?feddy"  area  on  the  downstream  surface.  Tem- 
perature differences  measured  by  Ranz  and 
Marshall  (1952)  show  a  similar  pattern  around 
a  ventilated  drop. 


Both  Frossling  (1938)  and  Kinzer  and  Gunn 
(1951)  introduce  "wind  factors"  by  which  the 
righthand  sides  of  [1]  and  [2]  are  multiplied 
to  account  for  increased  sublimation  due  to 
ventilation.  These  ventilation  factors  are  re- 
lated by  boundary  layer  theory  to  the  particle 
Nusselt  number  (Nu)  =  2rK/Dfor  heat  transfer, 
and  to  an  analogous  ratio  for  mass  transfer, 
the  Sherwood  number  (Sh),  often  denoted  by 
(Nu')  in  the  literature.  Thus,  the  theoretical 
equations  are 


(dm/dt)  =  2^Dr  (p  -p  )  (Sh) 


L  (dm/dt) 

s 


2-rrKr(Tr  -T)  (Nu) 


[31 


[4] 


where  (Sh)  and  (Nu)  must  be  determined  by 
experiment.  Theory  predicts,  however,  that 
(Nu)  is  a  function  of  Reynolds  number 
(Re)  =  2rv/v  and  Prandtl number  (Pr)  =  Cpy/K, 
while  the  analogous  Sherwood  number  is  a 
function  of  (Re)  and  the  Schmidt  number 
(Sc)  =  u/paD.  In  these dimensionless groups, 
V  is  the  ventilation  velocity,  fx  is  the  dynamic 
viscosity,  pa  is  the  density,  and  cp  is  specific 
heat,  all  for  air  in  this  case.  The  kinematic 
viscosity  is  v  =  y/pa. 

At  this  point,  the  Clausius-Clapeyron  equa- 
tion (i/Ps)  (dPs/dT)  =  lsM/rt2  may  be  used 
to  approximate  the  water  vapor  density  differ- 
ence in  [3],  where  Ps  is  the  saturation  vapor 
pressure,  M  denotes  the  molecular  weight  of 
water,  and  R  is  the  universal  gas  constant. 
Thorpe  and  Mason  (1966)  show  that  [3]  and 
[4]  may  be  combined  to  give  the  following 
expression  for  sublimation  rate  of  a  spherical 
ice  particle: 


dm 
dt 


2Trr  (p/P  -1) 
si 


1 


fL  M 
s 


KT  (Nu) 


RT 


-  1 


+ 


Dp  (Sh) 
s  I 


[51 


where  pST  represents  the  saturation  density  of 
water  vapor  at  temperature  T.  The  quantity 
(p/psT  -l)  is  the  undersaturation  of  the  en- 
vironment with  respect  to  water  vapor,  and. 
following  their  notation,   is   denoted  by  o  . 

Results  of  (heir  experimental  work  with 
ice  spheres  (0.3  mm  <  r  <1.8  mm)  at  ventilation 
speeds  from  25  to  100  cm  sec _1  show  the 
Nusselt  number  and  the  Sherwood  number  are 
practically  the  same,  and  are  proportional  to 
the  square  root  of  the  Reynolds  number. 

l  / 

(Nu)   =  (Sh)  =  1.88  +  0.58  (Re)  '2  [6] 

where  (Re)  ranged  between  10  and  200.  Values 


.2 


of  D  in  their  computations  are  about  7  percent 
lower  than  the  International  Critical  Tables 
value,  in  keeping  with  the  experience  of  other 
workers,  for  example,  Ranz  and  Marshall  (1952). 

Mean  ventilation  velocity,  V,  is  the  average 
relative  velocity  of  air  with  respect  to  the 
particle.  In  these  experiments,  particle  position 
was  fixed  in  a  moving  air  stream  of  low  tur- 
bulence intensity  (1  to  3  percent).  The  wind 
factor,  F,  by  which  sublimation  is  increased 
due  to  convection,  has  a  theoretical  value  of 
unity  when  the  Reynolds  number  is  zero.  How- 
ever, [6]  provided  the  best  fit  to  their  experi- 
mental data,  so  that  F  =  1/2  (Nu)  is  slightly 
less  than  unity  with  no  ventilation.  Ventilation 
may  increase  the  diffusion  of  water  vapor  from 
the  ice  sphere  by  as  much  as  five  times  if 
particle  Reynolds  numbers  as  high  as  200  are 
involved  (fig.  1). 
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Figure  1. — Experiments  by  Thorpe  and  Mason 
(1966)  show  that  ventilation  increases  sub- 
limation by  a  factor  (F)  proportional  to  the 
square  root  of  the  Reynolds  number,  where 
(Re)  =  (2rV/v)  ranged  from  10  to  200. 


3.  Particle  Temperature 

For  water  drops  evaporating  in  still  air, 
measurements  by  Johnson  (1950)  confirm  the 
relation  (Tr  -t)  =  (dls/k)  (p  -  pr)  .  This 
expression  for  the  temperature  difference  be- 
tween the  particle  and  its  environment  comes 
from  assuming  steady  state  equilibrium  of  heat 
and  mass  transfer  expressed  by  [1]  and  [2]. 
Multiplying  [1]  by  Ls  gives  an  expression 
equivalent  to  [2].  In  practice,  the  value  of 
p  r  is  taken  as  the  saturation  vapor  density  at 
particle  surface  temperature  Tr  .  The  value 
of  p  ,  the  ambient  vapor  density,  is  determined 
from  the  relative  humidity  of  the  environment 
and  its  temperature  T. 

Direct  measurements  by  Kinzer  and  Gunn 
(1951)  and  Ranz  and  Marshall  (1952)  show  that, 
for  practical  purposes,  the  surface  temperature 
of  a  water  drop  evaporating  with  ventilation 
is  equal  to  the  wet  bulb  temperature  estimated 
from  psychrometric  tables.  Kinzer  and  Gunn 
(1951)  also  found  the  temperature  of  suddenly 
ventilated  drops  adjusted  to  the  wet  bulb  tem- 
perature in  only  a  few  seconds. 

For  this  model,  particle  temperatures  are 
assumed  to  equal  ice  bulb  temperatures.  In  the 
temperature  range  from  -20°  C.  to  0  °  C.,  (Tr  -t)  , 
taken  as  the  ice  bulb  depression,  is  less  than 
0.5°  C.  for  relative  humidities  greater  than  70 
percent  over  the  elevation  range  0  to  4  km. 
Thorpe  and  Mason  (1966)  assumed  the  tem- 
perature difference  (Tr  -t)  was  so  small  for 
conditions  below  freezing  that  the  saturation 
vapor  density  at  the  particle  surface  was  taken 
as  that  for  ambient  temperature  T.  Thus, 
particle  temperature  does  not  enter  directly 
in  [5]. 


4.  Calculations 

Together,  [5]  and  [6]  provide  sublimation 
rates  for  ice  spheres  under  conditions  of  mod- 
erate ventilation,  at  temperatures  below  freezing, 
without  considering  the  radiation  of  heat  to 
the  particle.  Some  calculations  will  provide 
a  graphical  picture  of  the  relation  between 
factors.  Before  calculations  can  be  made  using 
[5],  however,  a  number  of  "constants"  must 
be  assigned  values.  This  has  been  done  in  the 
tables  in  the  appendix. 

Figure  2  shows  the  sublimation  rate  as  a 
function  of  particle  diameter  for  several  ambient 
temperatures  and  mean  ventilation  velocities. 
All  values  were  computed  for  1000  mb  atmos- 
pheric pressure ,  and  90  percent  relative  humidity. 
The  sublimation  rate  at  -10°  C.  is  almost  twice 
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Dingle,  and  Radok  (1966)  reported  mean  par- 
ticle diameters  near  100  \x  in  the  Antarctic, 
and  Kojima  (1969)  presented  particle  size  distri- 
butions with  modes  in  the  20-40  fx  diameter 
class,  measured  at  Hokkaido  University,  Japan. 
Measurements  on  water  drops  by  Kinzer  and 
Gunn  (1951)  showed  a  rapid  decrease  in  mass 
transfer  for  diameters  below  150  fx  when  the 
relative  velocity  was  equal  to  drop  fall  velocity. 
They  point  out  that  fog  and  cloud-size  drops 
(around  25  fx  )  evaporate  at  about  the  rate 
predicted  for  a  drop  at  rest.  At  low  Reynolds 
numbers  (0  to  7),  the  mechanism  is  not  com- 
pletely known,  but  air  is  apparently  entrained 
by  viscous  forces  and  the  drop  tends  to 
evaporate  into  an  air  shell  around  it.  (Radius 
of  curvature  does  not  exert  an  important  in- 
fluence on  vapor  pressure,  for  particles  larger 
than  1  micron.) 

Perhaps  the  best  limit  to  use  in  applying 
equation  [5]  to  small  particles  is  (Re)  =  10, 
the  lower  limit  given  by  Thorpe  and  Mason 
(1966)  which  agrees  with  the  arguments  of 
Kinzer  and  Gunn  (1951).  A  100  fx  particle, 
ventilated  by  a  100  cm/sec  air  stream  at  -20°  C., 
(Re  =  8.6)  is  slightly  outside  this  limit.  With 
90  percent  relative  humidity  and  1000  mb  atmos- 
pheric pressure,  equation  [5]  predicts  the  sub- 
limation rate  will  be    1.5 -10~9  g/sec. 


Figure  2. — The  sublimation  rate  of  a  single 
ice  sphere  almost  doubles  for  each  10-degree 
increase  in  ambient  temperature  between 
-20°  C.  and  0°  C.     Ventilation  rate   (V  in 
cm/sec)  has  a  smaller  effect,  according  to 
equation  [5]  with  values  at  1000  mb.  and  90 
percent  vapor  density  with  respect  to  satura- 
tion over  ice.      (The  dashed  portion  of  the 
curves  are  extrapolations  below  (Re)  =  10.) 


the  rate  at  -20°  C.  throughout  the  size  range, 
and  more  than  doubles  when  sphere  diameter 
is  doubled.  Examination  of  [5]  reveals  that 
sublimation  is  directly  proportional  to  the  water 
vapor  deficit  in  the  surrounding  environment. 
Doubling  this  undersaturation  doubles  the  sub- 
limation rate;  80  percent  relative  humidity  would 
give  twice  the  sublimation  rates  shown  in 
figure  2. 

Rate  of  sublimation  decreases  for  smaller 
particles,  but  the  loss  becomes  a  larger  per- 
centage of  particle  mass.  This  follows  from 
equations  [5]  and  [6],  which  predict  that  sub- 
limation rate  varies  as  the  3/2  power  of  par- 
ticle radius,  and  the  fact  that  total  particle 
mass  varies  with  the  third  power  of  radius 
(fig.  3). 

Particles  studied  by  Thorpe  and  Mason 
were  larger  than  most  windblown  snow.  Budd, 
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Figure  3. — The  sublimation  rate  calculated  by 
equation  [5]  represents  a    much  larger 
percentage  of  the  total  mass    for  small 
particles.     Calculations  plotted  here  are 
for  1000  mb  pressure  and  a  vapor  density 
that  is  90  percent  of  the  saturation  value 
over  ice.     Particle  density  is  assumed  to  be 
0.92  g/cm3. 
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5.  Radiation 

For  laboratory  studies,  most  experimenters 
have  justified  neglecting  heat  transfer  by  radia- 
tion as  part  of  the  sublimation  or  evaporation 
process.  However,  Bergen  3  shows  this  factor 
is  important  for  sublimation  of  a  snowflake 
under  solar  radiation.  Since  redistribution  of 
snow  by  wind  may  occur  with  clear  skies, 
radiation  will  be  included  in  this  model. 
Following  Bergen's  procedure,  a  term  is  added 
to  the  heat  transfer  relation  [2]  so  that 

L  (dm/dt)  =  4uKr(Tr  -T)  +  Q  [7] 

where  Q  is  the  net  rate  at  which  heat  is  trans- 
ferred to  the  particle  by  radiation.  If  the  heat 
transfer  expression  [7],  with  the  radiation  term, 
is  used  in  place  of  [2]  to  derive  the  equation 
for  particle  sublimation  rate  by  the  argument 
of  Thorpe  and  Mason   (1966),  the  result  is 


dm 

2iTra(Nu) 

-  _a 

KT 

L  M 

s     -  1 
RT 

dt 

L 

s 

fL  M 
s 

1 

+ 

KT 

RT 

DPST 

[8] 


To  evaluate  Q,  some  expression  for  the  total 
radiation  balance  on  the  particle  must  be  formu- 
lated. This  will  certainly  be  a  widely  varying 
term  since  conditions  range  from  clear  nights 
to  clear  days,  with  all  cloudy  and  foggy  situ- 
ations between. 

According  to  definitions  outlined  by  van  de 
Hulst  (1957),  a  particle  causes  extinction  of 
incident  radiation  equal  to  the  amount  it  scatters 
and  absorbs.  Scattered  radiation  includes  that 
defracted,  refracted,  and  reflected  by  the  particle. 
Of  the  incident  radiation  flux,  SQ  ,  the  portion 
absorbed  by  the  particle  is  equal  to  cabs  sQ  , 
when  cabs  is  an  area  defined  as  the  absorbing 
cross  section  of  the  particle.  An  absorbing 
efficiency  may  also  be  defined  so  thatEabS=cabs/G 
where  G  is  the  particle's  geometric  cross  section 
(G  =  irr2  for  a  sphere).  If  ;  is  the  particle 
albedo,   Eabs  =  1  -  £  . 

The  total  radiation  flux  intercepted  by  a 
spherical  particle,  suspended  in  turbulent  air 
above  a  horizontal  snow  surface,  will  be  taken 
as  the  sum  of  direct  and  scattered  solar  radia- 
tion. If  S0  is  the  total  solar  radiation  flux 
from  above,  and  as0  is  that  portion  reflected 
by  the  horizontal  snow  surface,  then  the  total 
radiation  load  on  the  particle  is  s0  +  xs0  > 
where  the  horizontal  surface  albedo,  "k  ,  may 
be  different  from  the  particle  surface  albedo. 
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The  amount  of  this  incident  flux  absorbed  by 
the  particle  is  cabs(S0+AS0) ,  or  EabsG(s0+xs0) . 
It  is  equal  to  Q,  the  net  rate  of  radiant  heat 
transfer  to  the  particle,  and  with Eabs  *  l  -  z, 
and      G  =  irr2     ,  for  a  spherical  particle 

Q  =  (1  -  +  A)ur2SQ  [9] 

Values  of  albedo  for  a  horizontal  snow 
surface  range  from  75  percent  or  less  for  old 
snow,  to  almost  100  percent  for  new  snow 
surfaces  (Geiger  1965).  A  conservative  value 
ofX=  0.8  has  been  used  in  calculations  below. 
No  estimates  of  snow  particle  albedo  could  be 
found  in  the  literature,  and  a  calculation  based 
on  particle  optics  was  not  formulated  because 
no  assumed  surface  characteristics  seemed  well 
founded.  Instead,  the  value  £  =  0.5  was  taken 
as  a  very  rough  estimate. 

This  analysis  neglects  thermal  radiation 
from  the  horizontal  surface  to  the  particle, 
because  large  differences  between  surface  and 
particle  temperatures  do  not  seem  likely  under 
the  natural  conditions  considered  here.  When 
snow  is  transported  by  wind  just  above  asphalt, 
or  a  bare  soil  surface,  thermal  radiation  might 
become  an  important  additional  heat  flux  to 
the  particle. 

Figure  4  compares  sublimation  rates  com- 
puted from  [8]  and  [9]  with  rates  predicted 
by  [5].  Most  obvious  is  the  fact  that  radiation 
contributes  more  in  sublimation  from  large 
particles.  The  sublimation  rate  by  [8]  and 
[9],  with  S0  =1.0  cal/cm2  min,  is  18  percent 
larger  than  when  radiation  is  neglected  for  a 
100  [x  sphere,  and  84  percent  larger  for  a  1  mm 
diameter  sphere. 


6.  Atmospheric  Pressure 

Rates  of  sublimation  at  different  elevations 
may  be  compared  by  evaluating  [5]  and  [8]  for 
the  corresponding  ambient  pressures.  Only  two 
parameters  in  equation  [5]  vary  significantly 
with  pressure.  Ventilation  factors  change  due 
to  the  variation  of  the  kinematic  viscosity,  and 
diffusivity  of  water  vapor  increases  with  de- 
creasing pressure. 

Variations  in  kinematic  viscosity  were  esti- 
mated from  the  defining  relation  v  =  p/pa 
where  fi,  the  absolute  viscosity  of  air,  is  almost 
independent  of  pressure.  The  density  of  the 
air  was  taken  aspa  =  p/cvRTv where  p  is  the 
pressure,  Cv  is  the  coefficient  of  compressi- 
bility, R  is  the  gas  constant  and  Tv  is  the 
virtual  temperature  (List  1968).  The  results 
(included  in  the  appendix)  show  that  Reynolds 
number,  and  thus  the  ventilation  factor  [6], de- 
crease with  decreasing  pressure.  ^The  variation 
in diffusivity'with pressure  is  D/D0=(T/T0)n(p0/p), 
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Figure  4. — Solar  radiation  is  more  effective  in 
increasing  the  sublimation  rates  of  larger 
particles,  according  to  this  evaluation  of 
equations   [8]   and  [9]   for  -20°  C.  at  1000  mb. 
The  vapor  density  was  90  percent  of  the  satu- 
ration value  over  ice,  and  a  ventilation  rate 
of  100  cm/sec  was  assumed  for  all  sizes.  The 
particle  surface  albedo  was  taken  as  X,  =  0.5 
and  the  horizontal  surface  albedo  as  \  =  0.8 
in  equation   [9] . 


where  D0  is  0.180  cm2  sec  _1  atpG  =  1000  mb 
(List  1968).  Computations  in  the  appendix  are 
for  a  temperature  of -20°  C. 

The  change  in  the  computed  sublimation 
rate  over  the  elevation  range  from  0  to  4  km 
(0  to  13,000  feet)  is  compared  for  two  particle 
sizes  in  figure  5.  Substantial  increases  in  sub- 
limation rate  are  predicted  for  higher  elevations, 
and  the  increase  is  more  significant  for  smaller 
particles.  The  increase  in  diffusivity  outweighs 
the  reduction  in  ventilation  rate  due  to  de- 
creased density  of  air  at  higher  elevations. 

Besides  changes  in  ventilation  factor  and 
diffusivity,  comparisons  employing  equation  [8] 
must  account  for  the  increase  in  direct  and 
diffuse  solar  radiation  at  higher  elevations. 
The  mean  radiation  flux  increases  by  two  to 
three  times  between  sea  level  and  an  elevation 
of  4  km  (Geiger  1965).  Figure  6  compares 
sublimation  rates  at  sea  level  (1000  mb)  with 
those  at  2.5  km  (750  mb)  and  4  km  (600  mb). 
Estimates  of  S  0  ,  the  radiation  flux,  are  listed 
with  the  corresponding  atmospheric  pressure. 


Of  course,  holding  temperature  and  humidity 
constant  over  this  elevation  range  is  not  realis- 
tic, but  figure  6  does  show  that  the  net  effect 
of  decreasing  atmospheric  pressure  is  an  in- 
crease in  sublimation  rate  at  higher  elevations, 
if  ambient  temperature  and  humidity  are  the 
same. 


7.  Turbulence 

Atmospheric  turbulence  levels,  under  con- 
ditions of  interest  here,  are  about  10  percent 
for  streamwise  velocity  fluctuations,  using  ap- 
proximations suggested  by  Lumley  and 
Panofsky  (1964).  (The  turbulence  level  is  the 
ratio  of  root-mean-square  fluctuating  velocity 
to  free  stream  velocity.)  Measurements  sum- 
marized by  [6]  are  for  turbulence  levels  of  1 
to  3  percent,  as  noted  before.  The  question  is 
whether  the  higher  natural  turbulence  level 
would  cause  a  greater  particle  sublimation  rate 
by  providing  more  ventilation  than  predicted 
by  [6]. 

If  a  snow  particle  followed  all  turbulent 
accelerations  of  air  motion  exactly  during  wind 
transport,  then  the  relative  velocity  would  just 
equal  the  particle  fall  velocity  due  to  gravita- 
tional acceleration.  However,  the  particle  to 
air  density  ratio  is  close  to  1000,  assuming 
the  windblown  snow  particle  has  a  density 
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Figure  5. — The  sublimation    rate  for  an  ice 
sphere  is     more  rapid    at  higher  elevations 
according  to  equation  [5]  evaluated  for  -20° 
C.  ambient  temperature  and  a  water  vapor 
density  equal  to  90  percent  of    the  satura- 
tion value  over  ice.     A  100  cm/sec  ventila- 
tion rate  was  assumed  for  both  particle 
di ameters . 
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Figure  6. --When     the     increase    of  solar 
radiation  at  higher  elevations  is  included , 
equations  [8]  and  [9]  provide  this  picture 
of  particle  sublimation  rates.  Ambient 
conditions  again  were  -20°  C.  and  90  percent 
of  water  vapor  density  at  saturation  over 
ice,  with  a  ventilation  rate  of  100  cm/sec 
for  all  particles . 

near  0.9  g/cm3  (ice).  Accelerations  required 
to  overcome  particle  inertia,  and  cause  it  to 
follow  the  air  velocity  fluctuations,  must  be 
created  by  additional  drag  forces.  These  can 
only  occur  if  the  relative  velocity  between  the 
particle  and  air  increases.  Thus,  the  instan- 
taneous relative  velocity,  uR  (defined  by 
uR  e  up  -  u,  the  difference  between  particle 
and  air  velocities)  is  a  fluctuating  vector  in 
turbulent  flow.  It  becomes  equal  to  the  par- 
ticle fall  velocity  (up  =  w)  when  the  air  is 
calm  (u  =  o). 

Friedlander  (1957)  showed  that  the  mean 
square  relative  velocity  characterized  convection 
effects  on  heat  and  mass  transfer  from  particles 
suspended  in  turbulent  flow.  By  making  some 
rather  restrictive  assumptions,  Tchen  (1947)  de- 
rived equations  of  motion  from  a  balance  of 
forces  on  the  suspended  particle.  Hinze  (1959) 
and  Soo  (1967)  have  both  summarized  his  work. 
A  generalized  solution  to  these  equations  of 
motion  (Chao  1964)  gives  the  mean  square 
relative  velocity  in  terms  of  the  mean  square 
fluid  velocity,  the  turbulent  energy  spectrum, 
and  a  quantity  determined  from  the  stochastic 
equation  of  relative  motion.  This  relation 
predicts  that  particles  of  a  given  size  will  have 
an  energy  spectrum  that  coincides  with  the 
energy  spectrum  of  fluid  motion  up  to  some 


frequency,  above  which  "slip,"  or  relative 
motion,  increases  with  increasing  frequency 
of  turbulent  fluctuations. 

With  the  energy  spectrum  suggested  by 
Lumley  and  Panofsky  (1964)  for  horizontal 
gustiness  in  the  atmosphere,  Chao's  procedure 
indicates  that  a  100  p  snow  particle  would 
completely  follow  all  turbulent  fluctuations. 
A  1  mm  particle  would  exhibit  "slip,"  that  is, 
it  would  not  follow  wind  fluctuations,  for  fre- 
quencies higher  than  1/10  cycle  per  second. 
Csanady  (1963)  derived  a  similar  equation  for 
mean  square  relative  velocity,  by  which  he 
calculated  that  a  0.5  mm  dust  particle  would 
follow  all  atmospheric  eddies  with  a  vertical 
dimension  greater  than  10  m. 

The  above  calculations  are  based  on  some 
rather  tentative  assumptions,  including  the 
absence  of  particle  interactions  and  particle 
spin  created  by  the  shear  flow.  More  restric- 
tive, drag  forces  are  assumed  to  obey  Stokes' 
Law,  which,  strictly  speaking,  would  limit  appli- 
cation to  situations  where  the  particle  Reynolds 
number  is  unity  or  less.  An  opportunity  and 
need  for  well  designed  experiments  under  atmos- 
pheric conditions  surely  exists  here. 

Turbulence  tends  to  increase  sublimation 
by  increasing  the  relative  velocity  above  the 
particle  fall  velocity.  However,  the  calculations 
described  above  do  not  indicate  that  this  in- 
crease is  large  for  the  normal  size  range  of 
windblown  snow  particles. 


B.  A  Nonspherical  Particle 

The  fine,  delicate  structures  called  snow- 
flakes,  created  during  precipitation,  are  certainly 
not  well  approximated  by  a  sphere.  As  noted 
in  the  introduction,  these  beautiful  forms  are 
quickly  destroyed  by  mechanical  forces,  when 
wind  begins  to  move  deposited  snow.  This 
action  reduces  both  the  fine  structure  and  the 
maximum  dimensions  of  the  particles,  so  that 
the  spherical  approximation  becomes  somewhat 
more  realistic.  The  photomicrographs  presented 
by  Budd  et  al.  (1966)  and  Kojima  (1969)  all 
show  windblown  snow  particles  without  any 
of  the  fine  dendritic  snowflake  structure.  How- 
ever, the  particles  have  irregular  surfaces  and 
shapes  and  are  more  often  oblong  that  spherical. 

By  the  analogy  between  water  vapor  diffu- 
sion and  electrostatics,  it  is  relatively  easy  to 
estimate  how  sublimation  rate  is  affected  by 
deviations  from  the  spherical  toward  more  ob- 
long particle  shapes.  The  calculations  in  this 
section  show  that  these  deviations  tend  to  in- 
crease sublimation  rates,  but  the  effect  is  not 
large. 
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1.  Diffusion  Shape  Factor 

According  to  McDonald  (1963),  the  first 
treatment  of  the  analogy  between  diffusion  and 
electrostatic  fields  should  be  attributed  to  Max- 
well, while  Houghton  was  first  to  use  the 
analogy  in  investigating  ice-crystal  growth. 
Equation  [1],  which  was  derived  from  this 
analogy,  can  be  modified  to  account  for  particle 
shape  by  replacing  the  radius  r  with  C,  the 
electrostatic  capacitance  of  a  conductor  with 
identical  size  and  shape.  Then  [1]  becomes 
dm/dt  =  AttCD(p  -  pr) .  McDonald  has  verified 
this  analogy  for  several  shapes  which 
also  have  a  theoretical  solution  in  terms  of 
electrostatics,  and  he  measured  the  value  of 
C  for  a  number  of  crystal  shapes  which  do  not 
have  theoretical  solutions. 

The  four  cases  in  which  C  may  be  deter- 
mined theoretically  are  (McDonald  1963): 

(a)  sphere,  radius  r,  c  =  r 

(b)  thin  circular  disk,  radius  r,  c  =  2r/ir 

(c)  prolate  spheroid,  major  and  minor  semi- 
axis  b  and  c,  c  =  A/in[(b  +  a)c]  where 
A  =  (b2  -  c2)h. 

(d)  oblate  spheroid  of  major  and  minor  semi- 
axis  b  and  c.  c  =  Ae/sin~le  where 
e  =  (1  -  c2/b2)%. 

The  limiting  value  for  a  thin  prolate  spheroid, 
c  =  b/in(2  b/c)  was  found  to  adequately 
represent  the  needle  crystal  form  where  b  is 
half  the  needle  length  and  c  is  the  radius  of 
the  midsection.  McDonald's  conclusion  was 
that  errors  from  approximating  the  true  value 
of  C  by  theoretical  values  are  less  than  many 
of  the  other  uncertainties  in  predicting  crystal 
growth  or  sublimation. 

The  sphere  has  the  smallest  surface  area  per 
unit  mass.  Any  deviation  from  the  spherical 
shape  will  increase  the  area-to-mass  ratio,  and 
will  result  in  a  larger  sublimation  rate  for  a 
given  particle  mass.  In  figure  7,  values  of  C/r 
are  plotted  as  a  function  of  the  ratio  of  major 
to  minor  semi-axis,  £  for  a  prolate  spheroid 
(formed  by  revolution  of  an  ellipse  about  its 
major  axis).  For  a  particle  with  the  same 
mass  as  a  sphere  of  radius  r,  the  rate  of  mass 
transfer  by  diffusion  would  increase  25  percent 
when  the  ratio  of  major  to  minor  semi-axis  is  5. 

2.  Sublimation  Equation 

When  we  consider  the  additional  factor  of 
particle  shape,  equation  [8]  may  be  modified 
to  provide  a  more  realistic  model.  Similar 
to  [3] ,  we  have 

dm/dt  =  4ttDCF     (p  -  p  )  [10] 


and  in  place  of  [7], 


dm         .  „__ 
L  - —    =  MCF^ 
s  dt  h 


(T    -  T)  +  Q 
r 


[11] 


where  Fm  and  are  the  ventilation  factors 
comparable  to  1/2  (Sh)  and  1/2  (Nu), respective- 
ly, and  the  radius  has  been  replaced  by  C,  the 
shape  factor.  Again  as  outlined  by  Thorpe 
and  Mason  (1966),  the  Clausius-Clapeyron  equa- 
tion is  used  to  approximate  the  vapor  density 
difference,  and  temperature  differences  are  as- 
sumed small.    The  result  analogous  to  [8]  is 


4ttC(P/psT  -  1)    _  (Q/KT) 


dm 
dt 


LsM 


RT 


-1 


s 
KT 


RT 


-1 


[12] 


+ 


si  m 


Ventilation  factors,  Fh  and  Fm  ,  must  be  evalu- 
ated by 
shape  C. 
these 


experiment  for  a  particular  particle 
From  the  work  of  Thorpe  and  Mason, 
factors  are  expected  to  have  the  forms 


Fh  - 


F  = 
m 


A  +  B   (Pr)     J (Re) 


l 


/■ 


A  +  B  (Sc)   ' d  (Re) 


[13] 


where  A  and  B  are  constants. 

The  radiation  term,  Q  cannot  be  used  in 
the  form  given  in  [9],  but  must  take  into  ac- 
count the  projected  area  corresponding  to  the 
particle  shape.  Thus  Q  takes  the  form 


Q  =  a  - 


+  (1  -OXA-S 
z  o 


[14] 


where  A  ±  and  A2  are  the  projected  areas  of 
the  particle  perpendicular  to  the  direct  and 
reflected  solar  radiation,  respectively.  Values 
for  A  ±  and  A  2  depend  on  the  average  particle 
orientation  with  respect  to  the  horizontal,  and 
on  the  sun  angle. 

3.  Comparison  with  a  Spherical  Particle 

Sublimation  rates  predicted  by  [8]  for  the 
spherical  particle  and  by  [  12]  for  the  nonspherical 
particle  cannot  be  compared  exactly,  because 
values  for  the  wind  factors  Fh  and  Fm  are, 
in  general,  not  known  for  nonspherical  par- 
ticles. However,  at  least  some  idea  of  the 
uncertainty  can  be  gained  by  considering  the 
extreme  case  of  an  infinite  circular  cylinder. 
Work  in  hot-wire  anemometry  (Hinze  1959) 
has  established  that,  for  a  circular  cylinder 
with  major  axis  at  right  angles  to  the  flow, 

1/  1/  1/ 

(Nu)  -  0.42(Pr)  /5  +  0.57(Pr)     3(Re)  '2 
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Figure  7. — The  sublimation 
rate  of  a  prolate  spheroid 
of  ice,  computed  from 
equation   [12],  is  plotted 
here  as  a  percentage  of 
the  rate  for  a  sphere  of 
equal  mass.     The  varia- 
tion of  C/r  (the  mass 
diffusion  shape  factor) 
is  also  shown. 


C  =  b/c 


At  temperatures  and  pressures  examined  here, 
the  values  of  (Pr)  are  near  unity,  so  that  the 
relation  for  an  infinite  cylinder  is  approximately 

(Nu)  a  0.42  +  0.57 (Re)  '2 

compared  to  the  value  from  [6]  for  the  spheri- 
cal case  (s  =  l) 

(Nu)  =  1.88  +  0.58 (Re)  12 

Assuming  the  oblate  spheroid  approximates  an 
infinite  circular  cylinder  when  the  ratio  of 
major  to  minor  axis  is  K  =  103  ,  the  value  of 
C/r  is  13.2  from  the  limiting  approximation 
c=b/ln(2b/c)  .  Then  if  the  value  of  (Nu)  varies 
as  C/r,  for  a  ratio  of  g  =  20  ,  the  Nusselt 
number  would  be  (Nu)  -  1.76  +  o.58(Re)^  , 
only  slightly  less  than  the  spherical  case.  The 
value  of  Reynolds  number  (Re)  is  taken  as 
2cV/v  ,  defined  by  the  minor  axis  of  the 
spheroid  (for  the  long  axis  perpendicular  to 
the  flow).  From  these  approximations,  it  ap- 
pears that  particle  shape  has  very  little  effect 
on  the  ventilation  factor. 

Again  taking  the  standard  conditions  for 
these  calculations  (T  =  -20°  C,  p  =  1000  mb, 
(RH)  =  90  percent)  with  a  steady  mean  venti- 
lation rate  of  100  cm/sec,  evaluation  of  [12] 
with  Q  =  0  provides  the  comparison  in  figure  7. 
For  a  ratio  of  major  to  minor  axis  of  5,  the 
sublimation  rate  without  radiation  is  only  about 
10  percent  larger  for  the  spheroid  than  for  a 
sphere  of  the  same  mass.  When  radiation  is 
added,  (q  +  0)  ,  the  effect  of  particle  shape 
on  sublimation  rate  depends  on  the  particle's 
orientation  with  respect  to  the  sun.  The  maxi- 
mum projected  area  of  a  prolate  spheroid  in- 
creases as  the  cube  root  of  eccentricity,  so  a 
spheroid  with    f  =  5,  oriented  with  its  long 


axis  perpendicular  to  the  sun's  rays,  inter- 
cepts 1.71  times  the  radiation  intercepted  by 
a  sphere  of  the  same  volume.  A  spheroid 
with  the  volume  of  a  100  u  sphere,  and  with 
{;  =  5,  would  have  almost  twice  the  sublima- 
tion rate  of  the  sphere,  when  oriented  so  that 
the  shape  effects  were  maximum.  However, 
if  the  particle  orientation  during  turbulent  sus- 
pension is  random,  which  seems  most  likely, 
the  shape  effect  on  radiation  interception  should 
disappear. 

C.  Summary 

The  objective  of  Part  I  was  to  explore 
the  relative  importance  of  numerous  factors 
that  determine  the  sublimation  rate  of  a  single 
particle.  Some  conclusions  drawn  from  the 
various  references  and  calculations  are: 

1.  Ice  bulb  temperature  is  a  reasonable  assump- 
tion for  the  temperature  of  a  sublimating 
snow  particle,  so  differences  in  particle  and 
air  temperatures  will  be  small  (less  than 
0.5°  C.)  under  natural  conditions  of  blowing 
snow. 

2.  The  sublimation  rate  nearly  doubles  for  each 
10°  C.  rise  in  ambient  temperature  in  the 
range  between -20°  C.  and  0°  C. 

3.  Particle  sublimation  more  than  doubles  when 
the  diameter  is  doubled.  The  percentage 
of  particle  mass  lost,  however,  is  greater 
for  smaller  sizes. 

4.  Sublimation  of  a  single  particle  is  directly 
proportional  to  undersaturation  of  the  en- 
vironment with  respect  to  water  vapor,  when 
radiation  is  neglected. 
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5.  Heat  transfer  by  solar  radiation  can  double 
the  particle  sublimation  rate  computed  for 
transfer  by  conduction  and  convection  alone. 

6.  The  increase  in  estimated  sublimation  rate 
at  higher  elevations  is  large  enough  to  be 
important.  For  the  same  temperature,  hu- 
midity, and  ventilation  velocity,  the  rate 
at  4  km  (13,000  feet)  is  about  30  percent 
more  than  that  computed  for  sea  level, 
neglecting  radiation.  When  solar  radiation 
and  its  increase  with  elevation  are  considered, 
the  rate  computed  for  4  km  is  approximately 
50  percent  higher  than  at  sea  level. 

7.  Atmospheric  turbulence  should  cause  sub- 
limation rates  higher  than  those  predicted 
from  laboratory  studies,  but  tentative  cal- 
culations indicate  this  effect  is  not  large, 
for  snow  particles  less  than  1  mm  in 
diameter.  Field  measurements  are  certainly 
needed  to  test  this  conclusion. 


8.  Particle  shape  appears  to  have  less  effect 
on  sublimation  than  most  of  the  factors  con- 
sidered. The  calculated  rate  for  a  spheroid 
with  a  ratio  of  major  to  minor  axis  as  large 
as  5,  is  only  10  percent  greater  than  for  a 
sphere. 

If  a  spherical  ice  particle,  ventilated  at 
100  cm/sec,  starts  sublimating  in  an  environ- 
ment with  90  percent  relative  humidity  at  1000 
mb  pressure  (sea  level)  and  -20°  C,  without 
any  radiation  transfer,  then  increases  in  radia- 
tion, elevation,  turbulence,  and  eccentricity  all 
tend  to  increase  the  sublimation  rate.  Even  if 
all  these  factors  are  neglected,  an  ice  sphere 
100  jj.  in  diameter,  under  the  conditions  above, 
has  a  computed  loss  of  mass  equal  to  approxi- 
mately 19  percent  of  its  total  weight  in  1 
minute.  This  calculation  indicates  that  the 
sublimation  of  blowing  snow  is  important 
enough  to  explore  further. 


PART  II 

SUBLIMATION  FROM  A  VOLUME  OF  AIR  AND  SNOW 


Having  considered  the  sublimation  of  a 
single  particle,  and  the  importance  of  factors 
that  govern  this  rate,  the  next  question  is, 
how  do  the  results  of  Part  I  apply  to  sub- 
limation rates  from  many  snow  particles  in  a 
volume  of  air?  The  two  major  problems  dis- 
cussed in  this  part  of  the  argument  are  (1) 
effects  of  particle  concentration  on  transfer 
of  heat,  mass,  and  momentum  within  the 
volume;  and  (2),  sublimation  rates  for  a  dis- 
tribution of  particle  sizes.  The  results  show 
that  even  maximum  natural  concentrations  are 
low  enough  that  particle  interactions  can  be 
ignored,  and  particle  sublimation  estimates  de- 
rived in  Part  I  apply  directly  to  each  particle 
within  the  volume  of  air  and  snow. 

A.  Concentration  Effects  on 
Particle  Sublimation 

The  expressions  in  Part  I,  for  sublimation 
of  a  single  particle,  were  based  on  assumptions 
that  humidity  and  temperature  gradients  ex- 
tended to  infinity.  Additional  particles  in  a 
unit  volume  present  sources  of  water  vapor 
closer  to  the  particle  in  question,  and  may 
thus  reduce  the  rate  of  mass  transfer.  As  con- 


centration increases,  the  drag  force  experienced 
by  the  particles  increases,  once  the  interparticle 
spacing  becomes  small  enough  that  the  fluid 
boundary  layers  around  the  spheres  interact. 
Surrounding  particles  represent  heat  sinks  that 
reduce  the  temperature  gradients  around  a  given 
particle,  and  radiation  scattered  by  other  par- 
ticles also  modifies  the  heat  transferred  to 
the  particle  in  question.  The  importance  of 
these  effects  on  sublimation  rate  depends  on 
the  concentration  of  particles. 

1.  Maximum  Natural  Concentrations 

Estimates  of  the  concentration  of  snow 
in  air  during  transport  are  available  from  the 
experiments  of  Budd  et  al.  (1966)  in  Antarctica. 
If  the  lowest  centimeter  of  snow  drift  is 
neglected,  concentrations  approach  the  mass 
density  of  air,  and  since  this  argument  is 
limited  to  evaporation  of  particles  in  suspended 
transport,  the  maximum  concentration  of  snow 
will  be  taken  as  M  -  pa  .  (Below  1  cm, 
transport  is  mainly  by  saltation  and  creep.) 
The  volume  fraction  of  particles  is  the  ratio 
of  volume  occupied  by  solids  to  the  total  unit 
volume,  denoted  by  (j>.  This  ratio,  multiplied 
by  particle  density  (  p    ),  gives  the  mass  con- 
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centration  of  snow,  M  =  pp  (j)  .  For  maximum 
conditions,  where  the  mass  of  snow  is  equal 
to  the  mass  of  air  in  a  unit  volume, (|>=pa/(pa+Pp). 
With  pa  =  l-icr3  g/ cm3  and  pp  =  0.9  g/cm3  , 
0  -  Pa/pp  >  so  only  about  one-thousandth 
of  a  unit  volume  is  occupied  by  snow  at  maxi- 
mum concentrations  above  the  one  centimeter 
level. 


2.   Correlations  of  Concentration  and  Transfer 
Processes 

Studies  which  correlate  mass  transfer  with 
voidage  (1  -  (j))  show  that  even  the  maximum 
concentration  of  snow  expected  above  1  centi- 
meter would  cause  no  important  reduction  in 
mass  transfer  due  to  particle  interaction.  These 
measurements  were  made  by  Chu  et  al.  (1953), 
Mullin  and  Treleaven  (1962),  and  are  sum- 
marized by  Soo  (1967).  Marshal  and  Langleben 
(1954)  expressed  the  effect  of  a  cloud  of  water 
drops  on  the  growth  of  an  ice  crystal  by  a 
term  equal  to  the  sum  of  all  drop  radii  in  a 
surrounding  unit  volume.  This  approach  also 
leads  to  the  conclusion  that  the  concentrations 
present  during  snow  transport  by  wind  do  not 
significantly  reduce  mass  transfer  from  that 
of  an  isolated  sphere. 

The  correlation  of  drag  coefficient  with 
Reynolds  number,  at  the  voidage  level  cor- 
responding to  maximum  snow  concentration, 
indicates  no  measurable  increase  in  drag  at 
Reynolds  numbers  less  than  200  (Soo  1967, 
fig.  5.2).  If  drag,  and  thus  relative  velocity, 
are  not  altered,  then  convection  effects  on  par- 
ticle sublimation  are  independent  of  snow  con- 
centration under  natural  conditions. 

JMo  studies  that  correlated  voidage  with 
conduction  and  convention  of  heat  were  found. 
However,  no  important  changes  in  these  modes 
of  heat  transfer  are  expected,  even  at  maximum 
concentrations  of  snow  in  turbulent  suspension. 
This  inference  depends  on  the  close  analogy 
between  mechanisms  of  heat  and  mass  transfer 
during  particle  sublimation.  The  nearly  identical 
values  of  Nusselt  and  Sherwood  numbers,  noted 
in  [6],  strongly  support  the  conclusion. 

Heat  transfer  by  radiation  is  a  completely 
different  process,  and  the  argument  reaches  no 
clear  conclusion  in  this  case.  If  each  snow 
particle  in  a  volume  of  air  receives  only  radi- 
ation from  outside  the  volume,  and  none  of  the 
radiation  scattered  by  surrounding  particles,  the 
process  is  called  "single  scattering,"  and  heat 
transfer  by  radiation  is  independent  of  particle 
concentration.  Multiple  scattering  occurs  within 
the  volume  when  the  particles  interact.  Based 
on  remarks  by  van  de  Hulst  (1957),  single 
scattering  probably  applies  over  most  of  the 


natural  snow  concentration  range,  and  can 
certainly  be  made  to  apply  in  calculations,  by 
limiting  the  dimensions  of  the  volume 
considered. 

In  summary,  concentrations  of  natural  wind- 
blown snow  particles  are  small  enough  that 
particle  interactions  can  be  ignored  for  estimates 
of  the  rate  at  which  snow  sublimates  in  a 
volume  of  air.  Thus,  within  the  limits  of 
optical  path  length  to  which  single  scattering 
applies,  the  sublimation  rate  for  a  volume  is 
simply  the  sum  of  individual  particle  rates. 
If  the  mass  of  snow  per  unit  volume  (M) 
is  composed  of  N  uniform  spheres,  each  with 
mass,  m,  sublimating  at  a  rate,  dm/dt,  then  a 
good  approximation  of  the  total  sublimation 
from  particles  in  the  unit  volume  would  be 
Ndm/dt  =  (M/m)dm/dt. 


B.  Effect  of  Nonuniform  Particle  Sizes 

Size  distributions  of  windblown  snow  par- 
ticles measured  by  Kojima  (1969),  Budd  et  al. 
(1966),  and  Lister  (1960),  show  that  sizes  range 
from  a  few  microns  up  to  300  n  or  more.  The 
more  recent  studies  also  show  skewness  of  the 
size  distributions  toward  smaller  particles.  Be- 
cause the  sublimation  rate  depends  strongly 
on  particle  size,  it  is  important  to  take  the 
distribution  of  sizes  into  account  in  this  model. 


1.   General  Relationships  for  Nonuniform  Sizes 

Suppose  a  unit  volume  of  air-snow  mixture 
contains  N  snow  particles  which  have  a  total 
mass  M.  A  nominal  particle  diameter,  x,  may 
be  defined  so  that  the  particle  mass  is  given 
by  m  =pp(ttx3/6).  Let  some  function  f  (x) 
describe  the  frequency  distribution  of  nominal 
diameters  in  the  unit  volume.  Then  the  num- 
ber of  particles  with  nominal  diameters  in  the 
range  x  to  (x  +  dx)  is  N/(x)dx,  the  total 
number  of  particles  times  the  relative  frequency 
for  that  increment  of  the  particle  size  distribu- 
tion. Multiplying  the  mass  of  a  particle  with 
diameter  x  and  the  number  of  particles  in  this 
size  increment  gives  the  mass  of  all  particles 
with  diameters  from  x  to  (x  +  dx) 
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Mx  =  ^p™  /6)N    ^x)dx  t15! 

Integrating  [15]  with  respect  to  particle  size, 
and  rearranging  terms,  the  total  number  of 
snow  particles  in  the  unit  volume  is 

N  =  M/[(p  tt/6)  f?    ^(x)x3dx]  [16] 
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This  relation  determines  the  number  of  particles 
in  a  unit  volume  of  air  and  snow  when  the 
mass  concentration  and  particle  size  distribution 
are  known. 

Let  (dm/dt)x  denote  the  sublimation  rate 
for  a  particle  with  nominal  diameter  in  the 
range  from  x  to  (x  +  dx).  Then  the  rate  of 
sublimation,  (dM/dt)  x  for  the  mass  associated 
with  this  particle  size  increment,  is  equal  to 
the  number  of  particles  of  this  size,  multiplied 
by  the  single  particle  sublimation  rate, 


(dM/dt) 


(dm/dt)  (x)dx 


[17] 


Summing  rates  from  all  particle  size  increments 
gives  the  sublimation  rate  for  the  unit  volume, 
so  that  by  integrating  [17]  with  respect  to 
particle  diameter, 


PC 

(dM/dt)  =  Nj0 


(dm/dt)     f (x)dx 


[18] 


This  is  a  general  expression  for  the  rate  at 
which  a  mass  of  nonuniform  particles  sublimates 
in  a  unit  volume  of  air  and  snow.  The  func- 
tional form  of  (dm/dt)  x  was  discussed  in 
Part  I.  Possible  forms  for  the  distribution 
function,  f  (x)  are  considered  below. 


2.  Specific  Size  Distributions 

Studies  that  report  size  distributions  for 
drifting  snow  are  too  few  to  assure  a  general 
distribution  function.  Several  reasonable  forms 
are  discussed  here,  and  evaluation  of  [18] 
with  the  two-parameter  gamma  function  sug- 
gested by  Budd  (1966)  provides  an  example. 

An  exponential  function  that  depends  on 
precipitation  rate  and  crystal  type  describes 
the  size  distributions  of  aggregate  snowflakes 
during  precipitation,  according  to  measurements 
by  Gunn  and  Marshall  (1958)  and  Ohtake  (1970). 
Within  a  few  hours  after  snowfall,  the  number 
of  windblown  snow  particles  in  a  size  class 
decreases  exponentially  with  increasing  size.  4 
This  work  shows  that,  after  several  hours  of 
snow  drifting  without  precipitation,  a  log- 
arithmic-normal function  provides  a  better 
description  of  the  drift  particle  size  distri- 
butions. 

Budd  (1966)  shows  that  the  log-normal 
function  gives  a  fair  approximation  to  the  size 
distributions  reported  in  Budd  et  al.  (1966). 


Schmidt,  R.  A.  Distributions  of  snow  particle  size 
during  wind  transport  over  a  ridge.  (Manuscript  in 
preparation  at  Rocky  Mt.  Forest  and  Range  Exp.  Stn., 
Fort  Collins,  Colo.) 


For  his  calculations  of  drift  content,  however, 
Budd  used  a  two-parameter  gamma  function 
which  gave  as  good  a  fit  to  the  size  distribu- 
tions, and  was  mathematically  more  convenient 
for  his  purpose.  Following  his  suggestion,  the 
distribution  of  particle  size  is  specified  as 


Kx)dx  =  3ar(a) 


-x/3  a-1  , 
e        x  dx 


[19] 


where  a  and  /3  are  the  parameters  of  the  dis- 
tribution, such  that  a/3  is  equal  to  the  mean, 
and  T{a)  is  the  complete  gamma  function. 
The  shape  parameter  is  a  ,  and  /S  is  the  scale 
parameter  (Thorn  1958). 

With  a  specified  particle  size  distribution, 
[16]  will  yield  an  expression  for  the  number 
of  nonuniform  particles  in  a  unit  volume.  Sub- 
stituting the  expression  [19]  for/(x)  in  [16] 
gives  the  number  of  particles  in  terms  of  mass 
concentration  and  particle  size: 


N  =  M/ 


(P  TT/6)f 
LP 


3ar(a) 


e-x/3x(a+3)-ldx 


[20] 


By  the  definition  of  a  gamma  variate  and  the 
fact  that  the  area  under  a  distribution  function 
must  equal  unity,  [20]  can  be  solved  to  give 
the  number  of  particles  per  unit  volume  in 
terms  of  the  mass  concentration  of  snow  and 
the  parameters  of  the  size  distribution, 


N  =  M(6/p  ir)  [r(a)/T(a  +  3)] 
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[21] 


The  sublimation  rate  of  this  particular  distribu- 
tion of  particles  is  determined  from  [18].  The 
integration  cannot  be  performed,  however,  until 
a  relationship  between  particle  size  and  venti- 
lation rate  is  formulated. 


3.  Particle  Size  and  Ventilation  Velocity 

A  ventilation  rate  equal  to  the  particle  fall 
velocity  was  assumed  in  Part  I,  when  the 
effects  of  turbulent  air  accelerations  were 
neglected.  Continuing  with  this  tentative  as- 
sumption, a  deterministic  relationship  between 
particle  size  and  fall  velocity  will  admit  inte- 
gration of  [  18] . 

Mellor  (1965)  summarizes  the  problem  of 
expressing  the  fall  velocities  of  drifting  snow 
in  terms  of  particle  size.  For  diameters  less 
than  60  /x  ,  particles  obey  Stokes'  Law  reason- 
ably well,  and  fall  velocity  is  proportional  to 
the  second  power  of  diameter.  In  the  range 
between  100  /xand  1  mm,  fall  velocity  is  directly 
proportional  to  the  first  power  of  diameter. 
Although  drifting  snow  particles  larger  than 
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1.5  mm  usually  are  not  transported  by  turbulent 
suspension,  such  particles  fall  at  speeds  pro- 
portional to  the  square  root  of  diameter. 
Recognizing  that  transitions  must  occur  between 
these  ranges,  Mellor  suggested  the  relation, 
w  =  w  0  exp(-xQ/x),  for  fall  velocities  of  drift 
particles  less  than  about  5  mm  diameter.  In 
this  expression,  wD  is  the  terminal  velocity  of 
large  particles  with  x  0  diameter. 

In  his  analysis,  Budd  (1966)  assumed  the 
fall  velocity  of  all  suspended  drift  particles  was 
w  =  C2x,  noting  that  this  relationship  would 
be  in  error  for  the  small  particles.  Depending 
on  particle  shape,  the  proportionality  constant 
C2  was  between  2440 /sec  and  3880/sec,  the 
latter  applying  to  spherical  particles.  This 
assumption  is  used  in  the  following  calculation, 
which  does  not  warrant  the  more  complicated 
relation  suggested  by  Mellor. 


4.     Calculation  of  Volume  Sublimation  Rate 


arctic  measurements  by  Budd  et  al.  (1966). 
The  mass  of  snow  (M)  in  the  given  volume 
would  be  0.01  g . 

If  the  distribution  of  sizes  is  approximated 
by  [19],  then  parameters  a  and  /3  must  be 
specified.  Budd  (1966)  obtained  values  of  the 
shape  parameter  (  a )  near  15.  Calculations 
here  are  for  several  mean  diameters(x  =  a$), 
with  a  =  15,  so  that  only  the  scale  parameter 
(/3)  varies.  These  distributions  are  plotted  in 
figure  8. 

Nighttime,  or  minimum  radiation  transfer 
(Q  =  0)  is  assumed,  and  particle  shape  effects 
are  neglected  so  that  equation  [5]  is  used  in 
place  of  [12]  to  express  single  particle  sub- 
limation. If  [5]  and  [6]  are  combined,  the 
expression  is 

^  =  C1[1.88x  +  0.58x  (Vx/v)1/2]  [22] 

where  the  mathematical  presentation  is  con- 
siderably reduced  by  the  substitution 


To  estimate  a  sublimation  rate  by  [18], 
first  the  total  mass  of  snow  in  the  volume  of 
air  is  required.  With  this  value  for  M,  and  a 
size  distribution  function,  such  as  [19],  the 
number  of  particles  in  the  volume  can  be 
calculated  from  [16].  Finally,  an  expression 
for  single-particle  sublimation,  such  as  [12] 
derived  in  Part  I,  is  combined  with  some  re- 
lation between  particle  size  and  ventilation 
velocity  to  give  an  integrable  expression  for 
(dm/dt)x  in  [18]. 

As  an  example,  suppose  the  volume  of 
interest  is  a  horizontal  layer  of  snow  and  air, 
1  cm  thick  and  1  m  square.  (The  reason  for 
considering  these  dimensions  becomes  apparent 
in  Part  III.)  If  this  layer  is  located  1  m 
above  a  horizontal  snow  surface,  then  a  snow 
concentration  of  1  g/m 3  might  be  expected 
for  15  m/sec  wind  velocities,  according  to  Ant- 
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Assuming  that  ventilation  velocity  is  equal  to 
particle  fall  velocity  (V  =  w),  and  further,  that 
fall  velocity  is  directly  proportional  to  particle 
diameter  (w  =  C2  x),  the  rate  of  sublimation 
for  a  spherical  particle  is 

i  /  j 

(dm/dt)    =  C . [1.88x  +  0.58  (C./v)     2x  ]  [24] 
x        1  I 

For  all  particles  in  the  unit  volume,  the  sub- 
limation rate  is  estimated  by 

(dM/dt)  =  [NC  /ear(a)]  /°°[1.88  x  +  0.58  [25] 
-L  o 

In     I    \l  1 2       2l     -X/&  a_1j 

(C^/v)     z  xz]e  x  dx 
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where  N  is  computed  from  [21]  and  C-,,  de- 
fined by  [23],  is  constant  for  a  given  set  of 
temperature,  pressure,  and  humidity  conditions. 

The  results  of  evaluating  [25]  for  the  several 
size  distributions,  sublimating  at  -20°  C,  1000 
mb  and  90  percent  relative  humidity,  are  shown 
in  figure  9.  A  value  of  C  2  =  3880/sec  was 
used  for  spheres,  and  particle  density  was 
assumed  equal  to  0.92  g/cm3  .  Sublimation 
rate  increases  markedly  when  the  same  mass  of 
snow  is  composed  of  particles  with  smaller  and 
smaller  mean  diameters.  Skewness  of  size  dis- 
tributions toward  smaller  diameters  has  the 
same  effect  of  increasing  surface  area  and 
thus  sublimation. 

An  average  particle  sublimation  rate  may 
be  computed  by  dividing  the  total  number  of 
particles  into  the  volume  sublimation  rate. 
Equating  this  value  and  the  righthand  side  of 
[24]  gives  the  diameter  of  a  particle  with  the 
average  sublimation  rate.  These  calculations 
were  carried  out  for  the  four  distributions  in 
figure  8  (table  1).  Average  sublimation  diam- 
eters are  nearly  the  same  as  the  distribution 
mean  diameters  for  these  examples. 

C.  Summary 

Concerning  the  effect  of  snow  concentration 
on  transfer  processes,  the  following  conclusion 
was  drawn: 

1.  For  heights  greater  than  1  cm  above  the 
snow  surface,  natural  windblown  snow  con- 
centrations are  too  low  to  have  an  effect 
on  heat,  mass,  and  momentum  transfer,  so 
that  the  relations  developed  for  a  single 
sphere  in  Part  I  apply  to  each  particle  in  a 
volume  of  drifting  snow.  It  follows  that  the 
total  sublimation  of  snow  in  a  volume  is 
the  sum  of  the  sublimation  rates  for  each 
particle  as  if  it  were  isolated. 
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Figure  9. — The  sublimation  of  a  given  mass  of 
snow  greatly  increases  as  the  mean  particle 
diameter  decreases  according  to  this  plot. 
Equation  [25]  was  evaluated  with  the  size 
distributions  in  figure  8  to  estimate  the 
sublimation  rate  of  0.01  gram    of  snow  in  a 
10^  cubic  centimeter  volume  at  -20°  C.  and 
1000  mb .     Water  vapor  density  was  again  90 
percent  of  saturation  over  ice. 


Introduction  of  a  particle  size  distribution 
produced  calculations  summarized  by  figure  9, 
which  led  to  the  conclusion  that: 
2.  Sublimation  of  a  mass  of  snow  in  a  volume 
of  air  is  extremely  sensitive  to  the  distribution 
of  particle  sizes,  due  to  the  surface  area 
relationship.    However,  there  is  apparently 
little  error  in  computing  sublimation  for  N 
spheres  with  the  same  diameter  as  the  mean 
of  the  distribution,  at  least  for  the  distribution 
function  and  sizes  considered. 


Table  1 . --Comparison  of  average  sublimation  diameters,  x 
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(x  in  y) 
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N 
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PART  HI 

SUBLIMATION  IN  A  VERTICAL  COLUMN  OF  DRIFTING  SNOW 


Natural  concentrations  of  drifting  snow  de- 
crease very  rapidly  with  height,  and  average 
particle  size  also  becomes  smaller  in  the  first 
few  meters  above  a  snow  surface.  When  these 
changes  are  expressed  as  a  function  of  the 
mean  wind  velocity  profile,  the  volume  sub- 
limation rate  (Part  II)  can  be  integrated  with 
height,  to  predict  the  sublimation  rate  in  a 
vertical  column  of  drifting  snow.  Evaluation 
of  this  general,  steady-state  equation  indicates 
that  water  losses  through  the  sublimation 
process  are  potentially  quite  important  in  a 
water  balance.  The  vertical  fluxes  of  heat  and 
water  vapor  required  to  maintain  a  steady  state 
provide  an  interesting  set  of  conditions  for  the 
vertical  profiles  of  temperature  and  humidity. 


A.  General  Equation  for 
Steady-State  Sublimation 

When  snow  is  drifting,  speeds  are  high 
enough  that  the  vertical  profiles  of  wind  in 
the  first  10  meters  above  a  horizontal  surface 
are  very  well  expressed  by 


u  =  (u^/k)   ln(z/z  ) 


[26] 


according  to  the  Antarctic  measurements.  The 
mean  horizontal  wind  component  (  u  )  is  pro- 
portional to  the  natural  logarithm  of  height 
(z)  measured  vertically  from  the  surface.  A 
friction  velocity^*  =  /x0/pa  ), defined  by  the 
ratio  of  surface  shear  stress  to  air  density,  and 
von  Karman's  constant  (k  =  0.4)  determine  the 
proportionality  factor.  The  integration  constant 
(z  0  )  is  often  called  a  "roughness  parameter" 
because  its  value  depends  on  conditions  of  the 
horizontal  surface. 

For  this  windspeed  profile  relation,  the 
eddy  viscosity  is  e  =  kulVz  .  Shiotani  and  Aral 
(1953),  and  Loewe  (1956),  assumed  the  eddy 
diffusivity  of  uniform  snow  particles  was  equal 
to  the  eddy  viscosity.  Then,  equating  the 
vertical  flux  of  snow  by  turbulent  exchange 
with  the  downward  flux  due  to  gravity  pro- 
vides an  expression  for  the  relative  number 
of  particles  at  any  height,  in  the  form 


N  /N     =  (z/Z)-w/kU* 
z  L 


particles  have  the  same  mass.  (Sommerfeld 
and  Businger  (1966)  presented  measurements 
which  cast  some  doubt  on  the  assumption 
that  particle  diffusivity  and  eddy  viscosity  are 
equal.  However,  as  Radok  (1968)  points  out, 
their  measured  particle  fall  velocities  seem 
unreasonably  high,  which  could  account  for 
the  discrepancy.  If  particle  "slip"  is  negligible, 
the  assumption  should  be  valid,  and  the  argu- 
ment will  continue  with  caution.) 

Budd  (1966)  assumed  that  particles  in  each 
small  size  class  were  distributed  with  height 
according  to  [27].  Photoelectric  particle  counter 
measurements  5  support  this  assumption.  When 
Nzx  denotes  the  number  of  particles  at  some 
height  z,  in  a  nominal  diameter  class  x  to 
(x  +   dx),  this  assumption  is  formalized  by 


N      =  N7  (z/Z) 
zx  Zx 


-w/kuj 


[28] 


Then  if  /(x)  is  the  distribution  function  of 
nominal  diameters  in  a  unit  volume  at  the 
reference  height  (Z),  the  number  of  particles 
in  a  size  class  at  that  height  is  NZx 
and  thus 


Nzf (x)dx 


N      =  Nfz/Z)    W/kU*  f(x)  dx 
zx  Z 


[29] 


For  all  particles  at  height  z,  with  nominal  diam- 
eters between  x  and  (x  +  dx),  the  sublimation 


rate  is  (dM/dt) 


N 


(dm/dt)  x  ,  or 


(dM/dt)       =  N   (dm/dt)    (z/Z)  w/kU*:P(x)   dx  [30] 
zx        Z  x 

Integrating  this  expression  with  respect  to  nomi- 
nal diameter  (x)  gives  a  relationship  for  volume 
sublimation  rate  at  height  z  as 

(dM/dt)     =  N         (dm/dt)    (z/Z)  W/KU*  f(x)dx[31] 
z        Z  J  x 

o 

The  final  step  toward  a  general  steady-state 
expression  for  sublimation  within  a  vertical 
column  of  drifting  snow  is  to  apply  [31]  to  a 
small  volume  with  vertical  increment  dz,  and 
integrate  with  height,  to  get  (dM/dt)  c  ,  the 
total  rate  for  the  column  as 


[27; 


(dM/dt)      =  N  JZ  r  (dm/dt) 

C  £    o       o  ^ 

(z/Z)~w/kU*f(x)dxdz 


[32] 


The  number  of  particles  (N  z  )  at  some  refer- 
ence height  (Z)  must  be  known  to  evaluate 
Nz  from  [27].  This  equation  also  holds  for 
the  mass  concentration  ratio  as  long  as  all 


5  Schmidt,  R.  A.  Distributions  of  snow  particle  size 
during  wind  transport  over  a  ridge.  (Manuscript  in 
preparation  at  Rocky  Mt.  Forest  and  Range  Exp.  Stn., 
Fort  Collins,  Colo.) 
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Evaluation  of  this  expression  requires  the 
following: 

a.  the  mass  or  number  concentration  at  the 
reference  height,  to  provide  Nz  ; 

b.  a  functional  expression  for  the  size  distri- 
bution, f(x); 

c.  at  least  two  points  on  the  wind  profile,  to 
determine  the  friction  velocity,  uA  ; 

d.  some  expression  for  particle  fall  velocity 
(w)  in  terms  of  the  nominal  particle  diam- 
eter (x); 

e.  an  integrable  relation  for  the  single  particle 
sublimation  rate. 

B.  Calculations  for  Constant 
Temperature  and  Humidity 

The  objective  of  evaluating  [32]  at  this 
point  is  to  demonstrate  the  potential  magnitude 
of  drift  snow  sublimation,  even  under  relatively 
cold,  moist  conditions.  Calculations  are  for  a 
vertical  column,  1  meter  square  in  horizontal 
cross  section  and  extending  from  1  cm  to  10  m 
above  the  horizontal  snow  surface.  Average 
windspeed  at  10  m  is  15  m/sec,  friction  velocity 
(u*)  is  52  cm/sec,  and  Z  =  1  m  is  specified 
for  a  reference  height.  As  in  Part  II,  the  mass 
of  snow  in  a  1  cm  layer  at  this  level  is 
assumed  to  be  0.01  g,  with  particle  sizes 
distributed  as  a  two-parameter  gamma  function 
where  a  =  15,  and  x  =  100  fi .  The  number 
of  particles  (N  z  )  in  this  layer  has  been  cal- 
culated as  17,220. 

Since  sublimation  depends  strongly  on  par- 
ticle surface  area,  and  therefore  particle  con- 
centration, the  vertical  drift  concentration  pro- 
file is  an  important  part  of  sublimation  cal- 
culations for  a  column.  A  general  equation  for 
the  mass  of  snow  in  a  layer  at  any  height  z 
may  be  derived  from  [29]  as 

M    =  N     (irp  /6)    r  (z/Z)"Wku*x3^(x)  dx  [33] 
z  Z  p  '0 

When  Budd's  assumptions  for  fall  velocity 
(w  =  C2  x)  and  size  distribution  [19]  are  made, 
the  drift  concentration  profile  is  given  by 

Mz  =  Mz  (1  -  C2g/ku^ln  z/Z)"(a+3)  [34] 

Figure  10  compares  the  average  of  profiles  meas- 
ured by  Budd  et  al.  (1966,  p.  95),  for  14.5 
m/sec  mean  windspeed,  with  the  concentration 
profile  computed  by  [34].  To  obtain  a  reason- 
able prediction,  the  value  C  2  =  2440/sec  had 
to  be  used.  The  value  C  2  =  3880/sec  (for 
spheres)  predicted  too  rapid  a  decrease  in  con- 
centration with  height.  (One  should  not  con- 
clude too  much  about  particle  shape  from  the 


comparison,  since  the  assumption  w  =  C  2  x 
only  roughly  approximates  the  true  relationship 
for  much  of  the  size  range.) 

An  estimate  of  the  mass  transfer  shape 
factor  (C  in  Part  I),  which  corresponds  to  the 
value  of  C  2  =  2440/sec,  can  be  derived  from 
drag  coefficient  -  Reynolds  number  correlations 
for  irregular  particle  shapes.  One  such  com- 
pilation of  measurements  is  presented  in  Report 
12  of  the  Subcommittee  on  Sedimentation  (Inter- 
Agency  Committee  on  Water  Resources  U.S.A., 
1958).  Comparisons  of  fall  velocities  predicted 
by  w  =  C  2  x  with  those  given  by  the  drag 
curve  for  different  particle  shapes  indicate  that 
the  capacitance  C  would  be  only  about  4  per- 
cent greater  than  for  a  sphere  (C  =  x/2),  even 
when  C  2  =  2440/sec.  The  conclusion  is:  even 
though  particle  shape  plays  an  important  role 
in  determining  the  drift  concentration  profile, 
no  large  error  in  sublimation  rate  is  introduced 
when  single  particle  rates  are  computed  for 
the  spherical  case. 

This  conclusion  leads  to  the  particle  sub- 
limation rate  equation  [8],  which,  with  [6], 
and  the  assumption  w  =  C  2  x,  again  results 


Mass  of  snow  in  l-cm  layers  of  a  column(g/IO  cm  ) 

Figure  10. — The    snow    concentration  profile 
computed  from    equation   [34]     with  = 
2440/sec  compares  reasonably  with  the  average 
drift  density  profile  reported  by  Budd  et 
al.    (1966)  for  windspeeds  near  14.5  m/sec. 
Values  used  in  the  calculations  included  a 
10-m  windspeed  of  15  m/sec,  u^  =  52  cm/sec, 
a  =  15  and  x  =  100  microns.     The  drift  den- 
sity at  the  reference  level  Z  —  1  m  was 
assumed  to  be  M    =  0.01  g/10^cm^. 
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in  [24]  for  particle  sublimation  without  radia- 
tion (Q  =  0).  When  this  expression  for 
(dm/dt)x,  and  the  distribution  function  [19], 
are  substituted  in  [31]  the  sublimation  within 
a  layer  at  height  z  is 


dM 
[dt 


=  Nzci 


z    ear(a)  J 


[1.88xa  +  0.58(C2/v) 1/2 


[35] 


x  a+1]eCC2a~1/e)xdx 


where  the  substitution  a  =  (1  ku*)  lnz/Z  has 
been  used  with  the  identity,  (z/z)_c2x/ku*  = 
exp  (c2ax) .    Evaluating  the  integral  gives 


dM 


dt 


=  NZCX  x(l  -  C2a6)   (a+1) [1.88  + 
0.58(C2/v)  /2(x+6)/(l-C2a6)] 


[36] 


The  final  step  of  integrating  [36]  with  respect 
to  height  does  not  lead  to  an  analytical  ex- 
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Sublimation  from  l-cm  layers  of  a  column  (g/sec) 

Figure  11. — Assuming  the  humidity  is  constant 
with  increasing  height  (in  this  case,  90 
percent  of  the  saturation  value  over  ice  at 
-20°  C),  the  sublimation  rates  in  1-centi- 
meter layers  of  a  column  of  air  and  snow 
show  a  variation  with  height  similar  to  the 
snow  concentration  profile.     Values  in  this 
figure  were  computed  by  equation  [36]  for 
the  calculated  concentration  profile  in 
figure  10. 
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pression  for  the  total  column  sublimation.  An 
approximation  has  therefore  been  obtained  by 
summing  the  sublimation  rates  computed  by 
[36]  for  each  of  the  layers  1  cm  thick  and  1  m 
square  that  make  up  the  column  of  drift  snow 
being  considered.  For  this  example,  the  varia- 
tion of  volume  sublimation  with  height  is  plotted 
in  figure  11,  which  illustrates  the  strong  influ- 
ence of  the  concentration  profile.  The  mass  of 
snow  computed  from  the  concentration  profile 
was  14  g  for  the  column,  and  the  estimated 
sublimation  rate  (dM/dt)c  was3.9*10~2  g/sec. 
Neglecting  the  strong  influence  of  radiation, 
even  at -20°  C.  and  90  percent  relative  humidity, 
the  estimated  snow  sublimation  rate  is  0.28 
percent  each  second,  which,  if  conditions  were 
somehow  maintained,  would  be  almost  17  per- 
cent per  minute. 


C.  Vertical  Fluxes  under  Steady  Conditions 

To  maintain  steady  humidity  and  tempera- 
ture conditions  within  the  column  of  drifting 
snow,  water  vapor  must  be  removed,  and  heat 
added,  at  rates  that  balance  the  sublimation 
process.  With  no  horizontal  gradients  of  vapor 
or  temperature  (steady  state),  turbulent  transfer 
of  vapor  to  the  atmosphere  above  the  column 
requires  a  vertical  humidity  gradient.  Heat  is 
added  to  the  column  from  the  snow  surface 
by  turbulent  exchange  along  a  vertical  tem- 
perature gradient  (forced  convection),  although 
both  direct  and  reflected  solar  radiation  may 
also  transfer  heat  to  particles  in  the  column. 
Thus,  radiation  and  the  vertical  temperature 
and  humidity  profiles  interact  with  the  sub- 
limation profile  to  establish  an  equilibrium  in 
this  steady  state  model. 


1.  Humidity  Profile 

If  the  humidity  of  some  thin  horizontal 
layer  in  a  column  of  drifting  snow  remains 
constant  with  time,  the  sublimation  rate  of 
that  layer  must  be  matched  by  an  increase 
in  the  vertical  water  vapor  flux  across  the 
layer.  Thus 


dM 
dt 


[37] 


where  E  denotes  the  vertical  vapor  flux  due  to 
turbulent  exchange.  Following  Priestly  (1959), 
the  vertical  flux  of  water  vapor,  E,  is  written 


E  =  -paKw(&r/8z) 


[38] 
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where  q  is  the  specific  humidity  (the  mass  of 
water  vapor  divided  by  the  mass  of  moist  air) 
and  K  w  is  the  eddy  diffusivity  for  water  vapor. 
Similarly,  for  momentum  flux  t  and  eddy 
viscosity  (Km)  the  defining  equation  is  x  =  pa 
Km(3u/3z) .  Specific  humidity  may  be  written  as 
q  =  p/pa  and  when  the  density  of  moist  air 
(pa)  is  considered  constant  over  the  column 
height ,  the  humidity  gradient  is(3p/3z)=pa(3q/3z). 
If  one  makes  the  basic  assumption  that  eddy 
diffusivity  is  equal  to  eddy  viscosity  (Kw  =  Km)  , 
the  vertical  vapor  flux,  in  terms  of  the  humid- 
ity and  windspeed  gradients,  is 


E  = 


p  (du/dz) 


Recalling  that 
be  rewritten 


E  =-uA(3p/8z)/OTT/3z) 


[39] 

,  the  flux  may 


[40] 


where  the  profile  of  specific  humidity  has  been 
written  in  terms  of  vapor  density  or  absolute 
humidity,  p  .  From  the  logarithmic  wind  pro- 
file, the  windspeed  gradient  is  (3u/3z)  =  u^/kz  , 
neglecting  the  zero  displacement  zQ.  So,  for 
conditions  of  interest  here  (near  neutral  stabil- 
ity), the  vertical  vapor  flux  is  estimated  from 


E  =  -u^kz  (3p/3z) 


[41] 


By  definition  a  e  (p/pgT-i) ,and  thus  3p/3z  =  psT 
3a/3z  so  [41]  may  be  written  as    e  =  -(ujcpgT)z 


(3a/3z)  and  differentiation  gives 


Jz~  =  -(u*kPsT> 


3^c 
3z< 


3a_ 
3z 


[42] 


The  equilibrium  condition  [37]  equates  [42] 
with  the  general  expression  [31]  for  (dM/dt)  z 
to  give  a  differential  equation  for  the  humidity 
profile. 


32a 


3a 


N. 


3z^      z  3z  ujcP 


sT' 


(dm/dt). 


[43] 


(z/Z)^/kU*  £(x)dx  =  0 


Single  particle  sublimation  (dm/dt)  x  also  con- 
tains the  variable  a  ,  and  for  the  steady  state, 
this  equation  is  an  ordinary  differential  equation 
of  second  order,  with  variable  coefficients.  If 
p,  Q_,  and  R  denote  functions  of  height,  and 
a  prime  indicates  differentiation  with  respect 
to  z,  [43]  has  the  general  form 

a"  +  ?o'  +  Qa    =  R  [44] 

where   R    =  0  when  radiation  is  neglected. 


Perhaps  because  of  the  unfortunate  combination 
of  assumptions,  no  solution  for  the  humidity 
profile  was  obtained  with  the  assumption 
w  =  C  2  x  and  the  two-parameter  gamma  dis- 
tribution for  f(x).  Nevertheless,  a  general 
picture  of  the  humidity  profile  curvature  is 
gained  from  the  vertical  vapor  flux  equation  [  41] . 

If  humidity  decreases  in  direct  proportion  to 
the  logarithm  of  height  (say  p  =  psT  -Alnz), 
then  the  humidity  gradient  is  inversely  pro- 
portional to  height  (3p/3z  =  -A/z)  ,  and  accord- 
ing to  [41],  the  vapor  flux  is  constant  with 
height.  However,  to  exhaust  the  vapor  created 
by  sublimation,  E  must  increase  with  height, 
to  match /z  (dM/dt) z  dz  as  shown  by  figure  12. 
Therefore"  the  humidity  profile  must  be  such 
that  z3p/3z  increases  with  height,  rapidly  at 
first.  Such  a  profile  is  sketched  in  figure  13, 
along  with  the  logarithmic  humidity  profile. 
Since  the  sources  of  water  vapor  are  distributed 
with  height,  the  vertical  humidity  profile  in 
drifting  snow  might  resemble  that  in  a  layer  of 
transpiring  vegetation. 
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Figure  12. — A  logarithmic  humidity  profile 
should  correspond  to  constant  vertical  vapor 
transport  with  height.     However,  the  total 
sublimation  increases  with  height  so  that 
such  a  humidity  profile  would  not  indicate 
the  required  balance  between  vertical  vapor 
flux  and  sublimation  rate.     The  humidity 
profile  used  in  these  calculations  was  az  = 
°100   (1  +  0.217 In  z/Z) ,  where  c^qq  is  the 
undersaturation  at  100  cm,   (taken  as  0.1). 
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Figure  13. — To    maintain  a  balance  between 
steady -state  vertical  vapor  transport  and 
sublimation  of  blowing  snow,  the  humidity 
profile  is  expected  to  be  similar  in  shape 
to  the  one  sketched  here . 


2.  Temperature  Profile 

For  the  vertical  heat  flux  (H),  a  defining 
equation  analogous  to  [38],  is  given  by  Priestley 
(1959)  as 


f- 


•p  c 

a  p 


_T 

3  Z 


[45] 


The  turbulent  heat-transfer  coefficient  is  Kh, 
and  the  adiabatic  lapse  rate  (  9  )  is  approxi- 
mately -0.01°  C/m.  Then,  with  the  logarithmic 
wind  profile  and  the  assumption  that  Kh  =K  m  , 
the  equation  becomes 


H  =  -p     C  u 
a    p  » 


kz 


+  G 


dZ 


[46] 


To  maintain  a  steady  sublimation  rate  within 
some  layer  at  height  z,  the  amount  of  heat 
required  is  Ls  (dM/dt)  z  .  If  this  heat  is  trans- 
ferred to  the  layer  entirely  by  turbulent  ex- 
change, the  vertical  heat  flux  is  reduced  by 
this  amount  across  the  layer  so  that  the  equi- 
librium condition  is 


M  =  -l  (dM/dt) 

dZ  S  Z 


47] 


equation  for  the  temperature  profile  in  a  form 
analogous  to  [43] . 


0+l    |T+  0 

3  Z  Z  3z 


p     C  u 
a    p  - 


_  1 
k  z 


dM 
dt 


=  0  [48] 


When  direct  and  reflected  solar  radiation 
are  included,  the  decrease  in  turbulent  heat 
flux  with  height  supplies  the  difference  be- 
tween heat  used  in  sublimation  and  that  trans- 
ferred to  the  layer  by  radiation  (Q  z).  Then 
the  equilibrium  condition  is 

3H/3z  =  -  [L     dM/dt  -  Qz]  [49] 

where  the  radiation  (Q  z  )  absorbed  by  particles 
in  the  volume  at  height  z  is 

Q     =  N„S     (1  +  A)  f  (1  -  r)   tt/4  x2  r(x)dx  [50] 

Z  Z    O  J0 

if  single  scattering  applies  throughout  the 
column.  The  temperature  profile  may  be 
obtained  from  the  differential  equation, 


3^T 
3^ 


z 


3z 


P-,      Cr,  U-^ 

a  p 


kz 


dM 
dt 


[51] 


P     C     u , 
a    p  ' 


kz 


Differentiating  [46]  provides  the  steady-state 


D.  Summary 

No  calculations  have  been  included  follow- 
ing the  arguments  on  the  humidity  and  tem- 
perature profiles.  Several  assumptions  should 
be  tested  experimentally  before  any  large 
amount  of  work  is  expended  in  such  calculations. 
Therefore,  the  assumptions  and  general  form  of 
the  model  are  listed  here  in  anticipation  of  the 
brief  discussion  that  concludes  this  paper. 
Important  assumptions  for  the  general  steady- 
state  equation  are: 

a.  The  sublimation  rate  of  each  windblown 
snow  particle  is  unaffected  by  particle  inter- 
actions. 

b.  Particles  of  each  size  class  are  distributed 
with  height  according  to  the  steady-state 
theory  for  uniform  particles  (where  eddy 
viscosity  is   equal  to  particle  diffusivity). 

c.  The  logarithmic  expression  [26]  describes 
the  vertical  wind  profile. 

The  sublimation  rate  within  a  unit  volume 
at  some  height  (z)  above  a  horizontal  snow 
surface  then  has  the  general  form  given  by 
[31]  as 

(dM/dt)     =  N     /"   (dm/dt)  (z/Z)~w/kU*/(x)dx 

Z  Z    J0  X 

A  differential  equation  for  the  vertical  profile 
of  humidity  results  from  the  assumptions  that: 

d.  All  sublimated  vapor  is  transferred  vertically 
by  turbulent  exchange. 

e.  Eddy  diffusivity  for  water  vapor  is  equal  to 
eddy  viscosity  (K  w  =  K  m ). 
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The  humidity  profile  is  obtained  by  the  solution 
of 

d2a      1  da  (dM/dt) 

dz^      z  dz   :  —     =  0 

uA  k  p  z 
sT 

Likewise,  for  temperature,  the  assumptions  are: 

f.  The  heat  required  for  the  sublimation  process 
is  supplied  by  radiation,  and  the  vertical 
transfer  by  turbulent  exchange. 

g.  The  coefficient  of  turbulent  heat  exchange 
is   equal  to  eddy  viscosity   (K  h  =  Km). 

The  differential  equation  for  the  vertical  profile 
of  temperature  is 

,2T       i    (AT         1     L   (dM/dt)  -Q 
d  T  +  1.    dT   |   Q      _s  z_    =   z 

dz2       z   Idz  p     C     u,  kz        p'    C    u.  kz 

'       a    p     ~  a    p  ~ 

DISCUSSION 

The  assumptions  used  to  derive  a  general 
relationship  for  sublimation  of  drifting  snow 
appear  quite  reasonable,  considering  the  low 
natural  concentrations  and  the  near-neutral 
stability  likely  under  these  strong  winds.  (The 
possible  exception,  in  assuming  that  eddy  vis- 
cosity equals  snow  diffusivity,  has  been  noted.) 
Many  of  the  suppositions  made  in  order  to 
evaluate  the  general  equation,  however,  present 
interesting  questions  for  the  experimenter. 

First,  what  is  the  density  of  a  drifting 
snow  particle?  Is  it  equal  to  that  of  ice,  as 
assumed  in  these  calculations?  If  not,  then 
the  surface  area  of  a  given  mass  of  drifting 
snow  is  larger,  and  the  particle  fall  velocity 
is  less,  than  predicted  here.  Answering  this 
question  is  fundamental  to  the  study  of  sub- 
limation from  drifting  snow. 

How  is  particle  fall  velocity  related  to  par- 
ticle size?  This  question  involves  the  answer 
to  the  particle  density  problem,  as  well  as  a 
correlation  between  particle  shape  and  fall 
velocity.  In  this  paper  the  linear  relation, 
w  =  C2  x,  was  used  to  express  the  fall  velocity 
in  terms  of  shape  and  size,  but  the  assumption 
has  an  important  influence  on  the  computed 
snow  concentration  profile.  Thus,  a  more 
accurate  approximation  would  be  worthwhile. 

What  is  the  actual  ventilation  velocity  of 
a  suspended  snow  particle?  Is  it  closely 
approximated  by  the  particle  fall  velocity,  or 
should  the  interaction  between  particle  size 
and  the  turbulent  energy  spectrum  be  included? 
Related  to  these  questions  is  the  effect  of 
ventilation  on  small  particles  (Re<10).  Cal- 
culations here  have  used  the  same  particle  sub- 
limation rate  throughout  the  size  range,  but 
this  may  overestimate  the  sublimation  of  the 
small  particles. 


By  the  arguments  presented  in  this  paper, 
the  vertical  profiles  of  both  humidity  and  tem- 
perature are  expected  to  show  the  influence 
of  sublimation.  Simultaneous  profile  measure- 
ments of  humidity,  temperature,  windspeed, 
and  snow  concentration  under  measured  con- 
ditions of  radiation  would  provide  a  test  of 
this  model.  However,  several  instrumentation 
problems  must  be  overcome  before  such  meas- 
urements are  possible. 

To  assess  the  hydrologic  importance  of  sub- 
limation from  drifting  snow  within  some  region, 
another  set  of  questions  requires  answers.  How 
frequently  do  strong  winds  move  snow?  What 
temperature  and  humidity  conditions  prevail 
during  these  events?  What  proportion  of  the 
drifting  occurs  at  night,  and,  during  the  day, 
how  much  solar  radiation  is  received? 

Probably  the  greatest  shortcoming  of  this 
sublimation  model  is  the  steady-state  assump- 
tions made  in  all  stages  of  the  development. 
Hopefully,  this  work  will  provide  some  initial 
level  of  understanding,  but  certainly  the  un- 
steady nature  of  the  phenomenon  must  be 
accounted  for,  in  a  true  description  of  sub- 
limation from  wind-transported  snow. 
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APPENDIX 


VALUES  OF  PHYSICAL  PROPERTIES  OF  AIR  AND  WATER  USED  FOR  CALCULATIONS 


Table  of  Values  at  1000  mb 


Temperature,  T  KJ_/  ^1/  ^(t)J/  v  yI/  M  r  M/r 


cal   cm      sec        C        ca 1  g        cm    sec  g  cm  cm    sec        g  cm      sec        g  mo  I         cal     K      mo  I        g     K  ca I 


-0° 

5 

80 

10"5 

677.0 

0.2Q6 

h 

8^7 

10"6 

0.1346 

1  .718 

10"'t     18.0153          1.98583  9.0721 

-5° 

5 

72 

10"5 

677.1 

.200 

3 

246 

10-6 

.1304 

1  .692 

10"* 

-10° 

5 

63 

ID'5 

677-5 

.194 

2 

139 

10"6 

.1259 

1 .667 

-15° 

5 

54 

10"5 

677-7 

.187 

1 

387 

10"6 

.1216 

1  .641 

lb"* 

-20° 

 5 

5^ 

ID"5 

677-9 

.180 

884 

10"6 

.1173 

1  .615 

ID"* 

Smithsonian  Meteorological  Tables  (List  1966.) 


International  Meteorological  Tables  (WMO-No.   188,  TP.  94);  taken  as  7  percent  less  than  table  value. 


Table  of 

Variation 

with  Pressure  at 

-20°C 

Elevation— 
(km) 

P 

2/ 

D— ^ 

S  5/ 
0 

mb 

kg/m3 

2  -1 

cm  sec 

2  -1 
cm  sec 

-2  . 
cm      m  1  n 

0 

1000 

1 .4039 

0. 1 150 

0.180 

0.8 

1 

900 

1 .2635 

.1279 

.200 

1  .1 

2 

800 

1 . 1231 

.1439 

.225 

1  .3 

750 

1 .0549 

.1532 

.240 

1.4 

3 

700 

.9867 

.1638 

.257 

1 .46 

650 

.9162 

.1763 

.277 

1.53 

4 

600 

.8457 

.1908 

.300 

1 .60 

—  Rough  estimate  from  NACA  standard  atmosphere. 
2/ 

—  Computed  from  virtual  temperature  (see  Smithsonian  Meteorological 
Tables  (List  1968) .) 


Computed  from  v  =  v/Pg- 

Computed  from  standard  equation  (see  text). 

A  roughly  parabolic  distribution  drawn  from  Geiger  (1965). 
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LIST   OF  SYMBOLS 


'abs 


C 

v 

D 

(dm/dt) 
(dm/dt) 

(dM/dt)c 

(dM/dt) 

(dM/dt) 

zx 

(dM/dt)  = 
z 

e 
E 

^abs 
*(x) 


a  =  -[1/kUj.  ln(z./Z)]  a  parameter  determined  by  the 

logarithmic  wind  profile  (sec/cm) 
b  =  major  semi-axis  of  ellipsoid  of  revolution  (cm) 
c  =  minor  semi-axis  of  ellipsoid  of  revolution  (cm) 
A  =  a  constant 

\^  =  project  area  of  particle  perpendicular  to  direct 
2 

solar  radiation  (cm  ) 
\?  =  projected  area  of  particle  perpendicular  to  reflec- 
ted solar  radiation  from  horizontal  snow  surface 
B  =  a  constant 

C  =  electrostatic  capacitance  (shape  factor  for  mass 
transfer) 

2 

particle  absorbing  cross  section  (cm  ) 
computational  grouping  of  factors  depending  on 

meteorological  conditions  (g/cm  sec) 

fall  velocity  shape  factor  (2440  sec      j<        jf^  3880 

sec  ~) 

specific  heat  of  air  [=0.240  cal/°K  (g  dry  air)] 
coefficient  of  compressibility 

2 

diffusivity  of  water  vapor  in  air  (cm  /sec) 
sublimation  rate  for  a  single  particle  (g/sec) 
sublimation  rate  for  a  particle  with  diameter  x 
(g/sec) 

total  sublimation  rate  for  an  arbitrary  column  of 
air  and  snow  (g/sec) 

sublimation  rate  of  mass  M  of  snow  in  an  arbitrary 
volume  (g/sec) 

sublimation  rate  of  the  mass  of  particles  M  with 
diameters  x  to  (x  +  dx)  in  an  arbitrary  volume  at 

height  z  Cg/sec) 
=  sublimation  rate  for  the  mass  of  snow  in  arbi- 
trary volume  at  height  z,  (g/sec) 
=  2.718  .  .  . 

=  vertical  flux  of  water  vapor  (g/cm  sec) 
=  particle  absorbing  efficiency 
=  general  particle  size  distribution 
=  1/2  (Nu) ,  wind  factor 
=  ventilation  factor  for  mass  transfer 
=  ventilation  factor  for  heat  transfer 

2 

=  particle  geometrical  cross  section  (cm  ) 


H 
k 

K 

\ 

K 
m 

K 
w 

L 

s 
m 

M 
M 


Zx 


Zx 

(Nu) 
P 

(Pr) 
P 
Q 

a 

r 
R 
R 

(Re) 
(BH) 


;  heat  flux  (cal/cm~sec) 

;  von  Karman's  constant,  0.4 

:  thermal  conductivity  of  air  (cal/cm  sec  °K) 
turbulent  heat  transfer  coefficient  (cm  /sec) 
2 

eddy  viscosity  (cm  /sec) 

2 

eddy  diffusivity  for  water  vapor  (cm  /sec) 
latent  heat  of  sublimation  (cal/g) 
particle  mass  (g) 

molecular  weight  of  water  (g/mol) 

mass  of  snow  in  an  arbitrary  volume  (g/volume) 

mass  of  snow  with  particle  diameters,  x  to  (x  +  dx) 

in  an  arbitrary  volume  (g/volume) 

mass  of  particles  with  diameter  x  to  (x  +  dx)  in  an 

arbitrary  volume  at  any  height  z  (g/volume) 

mass  of  particles  with  diameter  x  to  (x  +  dx)   in  an 

arbitrary  volume  at  reference  height  Z  (g/volume) 

exponent  in  diffusivity  equation  (1.75) 

number  of  particles  in  an  arbitrary  volume 

number  of  particles  in  an  arbitrary  volume  at  height  z 

number  of  particles  in  an  arbitrary  volume  at  the 

reference  height 

number  of  particles  with  diameter  x  to  (x  +  dx)  in  an 
arbitrary  volume  at  any  height  z 

number  of  particles  with  diameter  x  to  (x  +  dx)  in  an 

arbitrary  volume  at  reference  height  Z 

Nusselt  number  (2rK/D) 

atmospheric  pressure  (mb) 

Prandtl  number  (uC  /K) 
P 

general  function  of  height 

net  rate  of  heat  transfer  by  radiation  (cal/sec) 
general  function  of  height 
radius  (cm) 

general  function  of  height 
universal  gas  constant  (cal/mol  °K) 
Reynolds  number  (2rV/v)  >.r  (xw/v) 

relative  humidity  used  throughout  as  the  water  vapor 
density  in  percent  of  the  saturation  value  over  ice 


Schmidt  number  (u/p  D) 
a 


(Sc) 

(Sh) 

2 

S    =  total  solar  radiation  flux  (cal/cm  sec) 


Sherwood  number  (2rK/D) 
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T  = 


P 
UR 

~2 


time  (sec) 

environmental  temperature  (°K) 
particle  surface  temperature  (°K) 
temperature  of  horizontal  snow  surface  (°K) 
virtual  temperature  (°C) 
instantaneous  air  velocity  (cm/sec) 
mean  horizontal  wind  velocity  (cm/sec) 
instantaneous  wind  velocity 


(u  -  u)  relative  velocity 
friction  velocity  u 


o  a 


(cm/ sec) 


r(a) 


2  2 

mean  square  horizontal  velocity  of  air  (cm  /sec  ) 
ventilation  velocity  of  air  relative  to  particle 
(cm/ sec) 

particle  fall  velocity  (cm/sec) 

nominal  particle  diameter  [diameter  of  sphere  with 
equal  volume]  (cm) 
mean  particle  diameter 

height  measured  from  horizontal  surface  (cm) 
roughness  parameter  (cm) 
reference  height  (cm) 

parameter  of  particle  size  distribution 
parameter  of  particle  size  distribution  (cm) 
gamma  function  in  a 


.  e  =  eddy  viscosity  (cm  /sec) 
5  =  particle  surface  albedo 
X  =  albedo  of  horizontal  snow  surface 

y  =  micron,  10  ^  cm 

p  =  dynamic  viscosity  (g/cm  sec) 

2 

v  =  kinematic  viscosity  (cm  /sec) 

5  =  ratio  of  major  to  minor  semi-axis  (b/c)  for  prolate 
spheroid 

0  =  adiabatic  lapse  rate  (.01°C/m) 
TT  =  3.14   .    .  . 

3 

p  =  water  vapor  density  in  the  remote  environment  (g/cm  ) 
3 

o     =  density  of  air  (g/cm  ) 
a 

3 

Pp  =  particle  density  (g/cm  ) 

3 

p     =  water  vapor  density  near  particle  surface  (g/cm  ) 
psj  =  saturation  density  of  water  vapor  at  temperature  T 
(g/cm3) 

a  =  (p/pst  ~  1)  undersaturation  of  the  environment  with 

respect  to  water  vapor  [(RH)  -  100]  =  100a  in  percent 

2 

t  =  turbulent  shear  stress  (dynes/cm") 

2 

tq  =  surface  shear  stress  (dynes/cm  ) 
ij>  =  volume  fraction  of  snow  in  air 
-  <j>  =  voidage 

q  =  specific  humidity 
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